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Abstract 

We provide a Quantum Field Theory derivation of Lifshitz formula 
C3 ', for the Casimir force due to a fluctuating real scalar field in d + 1 

dimensions. The field is coupled to two imperfect, thick, plane mir- 
rors, which are modeled by background potentials localized on their 
positions. The derivation proceeds from the calculation of the vacuum 
^ ■ energy in the Euclidean version of the system, reducing the problem 

to the evaluation of a functional determinant. The latter is written, 
via Gelfand-Yaglom 's formula, in terms of functions depending on the 
structure of the potential describing each mirror; those functions en- 
code the properties which are relevant to the Casimir force and are the 
reflection coefficients evaluated at imaginary frequencies. 

1 Introduction 

Just a few years after Casimir's calculation of the attractive force between 
ideal metallic plates [1], Lifshitz derived his celebrated formula for the inter- 



action of bodies with arbitrary, frequency-dependent dielectric functions [2]. 
In the original calculation, Lifshitz considered two disjoint media-filled half- 
spaces with plane and parallel boundaries. The calculation was performed 
at finite temperature and the final result for the interaction force was writ- 
ten in terms of the dielectric functions that describe, macroscopically, the 
electromagnetic properties of each media. 

The impressive refinement achieved in precision experiments measuring 
the Casimir force stimulated a large amount of further theoretical calcula- 
tions. Indeed, to explain recent experimental results it has become increas- 
ingly important to use models that describe the mirrors in a more realistic 
way. For example, the dependence of the Casimir interaction on the geom- 
etry, temperature, and macroscopic electromagnetic properties of the inter- 
acting bodies is a subject of growing interest |3]. 

As an outcome of the sustained research effort on these topics, Lifshitz 
formula has been generalized in several directions. In this work, we shall focus 
on its generalization to the case of stratified media, i.e. a set of plane-parallel 
layers of different materials. The interaction between the different slabs can 
be computed and expressed in terms of the electromagnetic properties of 
the layers or, more specifically, in terms of their reflection coefficients, which 
generally depend on the field's frequency and momentum. 

Lifshitz formula, and some of its generalizations, have been derived using 
various different theoretical approaches. In the original setting [2], Lifshitz 
introduced a random field in the Maxwell equations, in order to describe 
the fiuctuating sources in the materials. The interaction between bodies was 
computed by evaluating the appropriate component of the stress tensor, after 
assuming a local correlation for the random field. Besides, the existence 
of different media was taken into account by imposing the corresponding 
boundary conditions for the Green's functions of the electromagnetic field. 
This formula has been recently rederived [3] using a canonical quantization 
approach for macroscopic QED, in which the starting point is the action of the 
electromagnetic field, coupled to a set of reservoir oscillators. An alternative 
approach |5], is to quantize the electromagnetic field in the presence of the 
plane-parallel layers, computing the zero-point energy by following Casimir's 
original line of thought. Then the boundary conditions on the modes of 
the electromagnetic field determine the allowed eigenfrequencies, and the 
summation of the zero-point energy of each mode is performed using the 
argument theorem. Note that, in this derivation, the validity of Lifshitz 
formula in lossy media is not apparent. 

Lifshitz formula has also been derived in a quite different context, namely, 
considering the mirrors, and the spaces between them, as a quantum optical 
network [6]. Within this formalism it becomes possible to compute the force 



between lossy mirrors, and the final answer can be expressed just in terms of 
the frequency-dependent reflection coefficients of each mirror. Remarkably, 
these coefficients are the only properties of the media that become relevant 
to the computation of the vacuum force, in spite of the fact that the same 
coefficient may be obtained from different kinds of mirrors. 
The Casimir force for absorbing media has also been considered in Ref. [7], 
where Lifshitz formula has been derived using the theory of quantum open 
systems. 

Another route to the computation of the Casimir force is to consider 
a vacuum field in the presence of background potentials localized on the 
mirrors [8], [9]. These models can be justified from a microscopic point of 
view taking into account the interaction of the internal degrees of freedom 
of the mirrors with the vacuum field P, [TU]. Using a functional approach, 
the integration of the internal degrees of freedom of the mirror produces 
an effective action for the vacuum field. That effective action contains a 
potential that is different from zero just at the positions of the slabs. These 
models may be used to reproduce different boundary conditions for the fields 
which approach perfect ones only under certain specific limits. Although 
localized on each mirror, the potentials are in general nonlocal in time as well 
as in the coordinates parallel to the mirror. In other words, the boundary 
conditions for the Fourier modes of the field may depend on frequency or 
momentum. 

In this kind of approach, the zero-point energy is given by the functional 
determinant of an operator containing the potentials that describe the me- 
dia, and the formal general expression seems to be strongly dependent on 
the specific shape of the potentials involved. This is however, a puzzling 
circumstance. As already mentioned, Lifshitz formula is a quite general ex- 
pression for the Casimir force between two plane and parallel mirrors, and 
the force depends on the potentials only through their reflection coefficients. 
In this paper we show, by providing a field theoretic derivation of the Lifshitz 
formula, that that is indeed the case. 

The derivation relies upon the use of Gelfand-Yaglom formula for functional 
determinants [TT] (for a modern review, see [12]). As a by-product, we shall 
find an expression which also holds true in situations where an uncritical ap- 
plication of Lifschitz formula may be problematic, namely, potentials leading 
to bound states. 

This paper is organized as follows: in Section |2] we define the kind of 
system we consider and derive an expression for the Casimir force. Then 
in Section [2] we discuss some consequences of the general result, considering 
some particular examples. Our conclusions are presented in Section |H 



2 Derivation of the Casimir force 

Throughout this paper, we shall consider the case of a single, massive, real 
scalar vacuum field, </?, in d + 1 spacetime dimensions, equipped with an 
Euclidean action S{ip), which has the structure: 

S{v) = Soiv) + Si{^) , (1) 

where 5*0 is the free field action 

5oM = lJd'^'x[{d^r + mV], (2) 

while Si describes the interaction between the field and two mirrors. Those 
mirrors are assumed to be plane and parallel, with their normals pointing 
along the x^ direction (the remaining directions will be denoted by Xq and 
x\\ = Xi, ...Xd-i). Their interaction with ip shall be described here by a 
potential V, local in x^, which, in view of the previous assumptions, will be 
a function of Xd concentrated on the regions occupied by the mirrors, and 
vanishing elsewhere. Thus the form of Si shall be: 

Siiip) = ^ j d'^+^xd'^+'x' 6{xd-x'd)V{xd,Xo-x'„x\^-x[)^{x)ip{x') , (3) 

which, as well as 5*0, is invariant under translations in all the spacetime 
coordinates except Xd, namely, under the shifts x^ -^ x^ + c^ with /i = 
0, 1, ... , Xd-i, c^ = constant. Note that we have included a non-local depen- 
dence of the potential on the temporal and parallel coordinates, in order to 
describe more general responses of the mirrors. 

The fact that V is concentrated around each mirror, for Xd G [oi, ^i] and 
Xd G ['^2; ^2]; say, may be made explicit by writing it in terms of two functions, 
Vi and V2, with support on the intervals [0,6i] and [0,^2], respectively, such 
that: 

V{xd,Xo-XQ,xii-x\^) = Vi{xd~ai)\i{xo-XQ,xii-x\^) 

+ V2{xd - a2)A2(xo - Xo,X|j - xj|) . (4) 

The functions Aj characterize the response of each mirror. We regard the 
distances 6i = bi — ai, i = 1, 2 as the 'sizes' (i.e., widths) of the mirrors, while 
I = a2 — bi is the distance between them. 

To give sense to the forthcoming steps, we confine the system to a d 
dimensional spatial box (containing the mirror) such that the field satisfies 
Dirichlet conditions on all the 2d boundaries. That box is assumed to have 



sides of equal length, Ly, for the d — 1 coordinates which are parallel to 

the plates: —-2^ <a;i < -^, i = 1, . . . ,d — 1, while for the remaining spatial 
coordinate, \xd\ < f, with L not necessarily equal to L||. Besides, the Xq 
coordinate is also assumed to have a finite range, Xq G [— T/2,T/2], and the 
field to satisfy periodic boundary conditions on that interval. 

We have represented, in Figure 1, our previous conventions and notations 
about the mirrors' configurations, from the point of view of the Xd coordinate. 
Regarding the shape of the potentials, the usual, potential barrier case has 
been depicted in Figure 1, since that it the case when considering imperfect 
Dirichlet conditions. However, the derivation below is independent of that 
assumption, and indeed, one can even consider potential wells. 




Xd 



Figure 1: The potential V{xd), using a typical profile for the mirrors. 

The vacuum energy per {d— l)-dimensional mirror volume, £, relative to 
the vacuum energy in the absence of the mirrors (V^ = 0), for a finite L, can 
be expressed in terms of two Euclidean vacuum persistence amplitudes: 



8 = - lim 

2 r,L|i->oo 
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(5) 



where Z, the vacuum amplitude corresponding to S, can be written as the 
functional integral: 

Z = fv^e-^^'^^ , (6) 



while Zq = Z\y^q. Both path integrals are performed over the space of fields 
vanishing on the boundaries of the spatial box. 

The exact -albeit formal- solution for the Z integral may then be given 
in terms of a functional determinant: 

Z = (detr)"^ , (7) 

with T = —d^ + m? + V , an operator acting on functions of x G W^^^ which 
vanish on all the boundaries of the box, and d^ = 8^8^, fi = 0,1, . . . ,d. 

Separation of variables allows one to reduce the problem to another one 
for an operator acting on functions of just one variable, x^, but dependent 
on a momentum vector k = (/cq, . . . , kd-i). The determinant becomes then 
block-diagonal: 

det r = W det f{k) (8) 



with 7'{k) = — 9J + Q,'^(k) + V{xii, k), VL{k) = y/k'^ + m^, which now acts on 
functions of Xd which vanish at x^ = ±-|. For a potential of the form given 
in EqSlwe have 

V{xd,k) = Vi{xd-ai)Xi{k) + V2{xd-a2)X2{k) = Vi{xd,k) + V2{xd,k) . (9) 

The values of k become continuous and unbounded for the limit in ([5j), 
which we take now: 

1 f d'^k ^ rdetf(A;)i , , 

S = - -rlog ^^^ , 10 

2 J {2nr ndetro(A;)J ^ ' 

with % = f\v=o- 

It is interesting, at this stage, to note that the only change one should 
make if one wanted to calculate the free energy per unit volume, F, for the 
same system at a finite temperature /3~^ would be to keep the xq coordinate 
finite, with T = fi. The integral over ko would then become a sum over the 
discrete Matsubara frequencies Un = %^. Thus 



If d'^^'-ki] 



detTioJn, k\\ 

V n, I!. q^ 



^Kt^J (2-)^-^ "LdetroK,fc||) 

where k\\ = {ki, . . . , kd^i). Nevertheless, because the only important differ- 
ence between the evaluation of the energy or the free energy emerges after 
the evaluation of the integrand, we just consider the energy, commenting on 
the analogue result for the free energy at the end. 



Note that the energy £, besides being defined for a compact Xd coordinate 
i\xd\ < ■§)) still contains a contribution from the mirrors' self energies. One 
possible way to get rid of those contributions is to consider the force (per 
unit volume) J^{L), and its L — )■ oo limit, the usual, L-independent Casimir 
force, J^c = J^{oo): 

d£ 

:f{l) = -7^ , :fc = lim :f{l) . (12) 

To find the ratio between determinants appearing in ([5j), we apply Gelfand- 
Yaglom's (GY) theorem, which (after a rescaling and shift of coordinates) 
allows we to write it as follows: 



det r{k) ^( 



•L^ 



(13) 



det7^(A;) ^0(2/ 
where ip and ipQ are solutions, respectively, of the homogeneous equations: 

fik)iP{xd) = , foik)M^d) = , (14) 

such that ijj{—^) = 0, V^'(— ■§) = 1, and identical conditions for ipQ. 
Coming back to (I12p . we note that. 



^L) = -- 7^^7^1og[:^], (15) 



1 f^^d_ jKf) 

2 J {2nYdl °^^^o(f)^ 

the Casimir force obtained afterwards by taking the L -> 00 limit. 

In fact, because of the logarithmic derivative above, any / dependent 
factor may be dropped, in particular t/^q, which is independent of V , and 
hence of /. 

Thus the next step is the calculation of '?/'(-|) from the second order dif- 
ferential equation; in what follows, we shall use x to denote the Xd variable, 
since the problem is essentially one-dimensional. We first note some proper- 
ties of the solutions to the homogeneous equation above, without making any 
assumption yet about V , except that it is regular enough as to make the ex- 
istence and uniqueness theorem for the second order homogeneous equation 
valid. Thus, given the values of ip and its derivative ip' at a single point, Xj, 
say, their values at any other point, x/, shall be uniquely determined. Intro- 

ducing the two-component x-dependent vector ^(x) = I /// wq ); where 

we divided by Vt in order for the components to have equal dimensions, after 
some (linear) algebra we may write: 

^(x/) = A{xf,Xi)^{xi) (16) 



where the matrix A can be written as follows: 

A{xf,x^) = U{xf)U-\xi) (17) 

with 

TT(^^ _ f Xi{x) X2{x) 

^^"^^ [x[{x)/n ^,{x)in 

Here, Xi ^"^^ X2 denote two independent solutions of the homogeneous equa- 
tion, forming a basis of solutions (it can be shown that A is independent 
of the choice of basis). From the fact that these two functions are indepen- 
dent solutions of the homogeneous equation, it follows that detU{x), their 
(scaled) Wronskian determinant, is a non vanishing constant. Therefore, we 
obtain the property: 

det A{xf,Xi) = 1 . (19) 

Note that ip{^) may be written in terms of a single matrix elements of 
A, as follows: 

V'(^) = Au{L/2,-L/2), (20) 

where we dropped / independent factors. 

For the particular case we are considering, namely, a potential which van- 
ishes everywhere except at the region occupied by the mirrors, the previous 
expression may be rendered in a more explicit fashion. To that order, it is 
convenient to make repeated use of (IT6|) to write A as a product of matri- 
ces corresponding to subintervals. Those factors are completely determined 
by the potential inside the respective intervals; thus, it is either just an Aq 
factor, where the potential vanishes, or it is determined by the solution of a 
homogeneous equation with V replaced by either Vi or V2. Thus: 

AiL/2,-L/2) = A^'\L/2,b2)A(-'''^\b,,a2) 

X A(°)(a2,6i)A(^^H6i,ai) AW(ai,-L/2) (21) 

where we have indicated which potential each factor corresponds to. The 
matrix A^^\xf,Xi), corresponding to a null potential on the interval [xj,X/] 
can be immediately found, by using, for example, the independent functions 
Xi{x) = e^^ and X2{x) = e~^^, the result being: 

.(».(.„., . ^<»),A.) ^ ( -^^ ^;^ ) , (22) 

Ax = Xf — Xi. Besides, note that the only nontrivial factors are A^^^'^'^ = 
^(i>2)^ 2x2 matrices depending on each potential, the wave vector k; they 
depend on a,, 6, only through their differences, 6i. 

8 



Taking into account the results above, we may write a quite compact 
expression for ip{L/2): 

ij{L/2) = u^A(2)a(°)(/)A«i;i, (23) 

where we introduced the 2-component vectors, Uj Vi, 

( cosh(fi/i) \ ( sinh(n/i) \ 

""' - \ sinh(n/i) J ' ""' - \ cosh(fiZi) ) ' ^^^> 

and /i, I2 are the distances from the mirrors to the boundaries of the spatial 
box. We shall assume now that hi = —1/2 and 02 = 1/2, so that the internal 
faces of the mirrors are symmetrically disposed with respect to the origin (as 
well as the boundaries at ±L/2). Then li^2 = {L — l)/2 — (5i^2- Note that, in 
(123|) . the only factors that depend on / are 1*2, fi, and A^^\ 

To calculate the logarithmic derivative of ijj{L/2), we first note that 



dij{L/2) ^ ^ 



dl 






2 



where 



r(o) ^ f smh{ni) cosh(fi/) \ , . 

V cosh(fi/) sinh(fi/) J ■ ^ ' 

Rather than writing the rather lengthy expression for the logarithmic deriva- 
tive for a finte L, we directly present its y — )■ 00 limit. Introducing the 

(constant) vector uii = -j^ [ 1 J > and matrix C = i _. ^ j : 

r ^log^P{L/2) ^ .^, wjA(')CA^'->wi 

L dl J^^- wfAi^^A^AWwi ■ ^ ' 

It becomes clear from the last expression that it is convenient to perform a 
change of basis in all the matrices involved. Indeed, rotating to the basis Wi, 

W2, with W2 = ^ \ 1 , and denoting by T^^'^'^^ and D the form adopted 

by A^'^'^''^^ and C in the new basis, respectively, we see that: 



More explicitly, 






dl 



rp{2) „(1) 

1 _i_ -'12 -^21 ^-2rz< 



o ^ rp{2) rp{l) 



-^12 -^21 -2ra 

7^(2)^.(1)^ 
-^11-^11 



Finally, we arrive to an expression for the force per unite volume. 



I f d'^k d 



-^=-27 (2^9/^°^ 



(2) m(l) 






12 -^21 -2™ 



-^11 -^11 



(29) 



(30) 



which is the main result of this paper. 

It is worth interpreting, at this point, the meaning of the ratios between 
matrix elements of T^^'^-* appearing in (130|) . From the original definition of A, 
and implementing the change of basis, one sees that the matrix elements of 
T(*\ 2 = 1,2 relate ?/; and il)' just to the right of the mirror with the function 
and its derivative just to the left, but in a mixed fashion. Indeed, in the new 
basis, it relates functions such that iIj' = Qtp or tp' = —Qtp from one face of 
the mirror to the other. Thus, the matrix elements of T^''' connect solutions 
of the form e^^^ on both sides of the mirror. 

We can make contact with the usual expression of Lifshitz formula from 
the following observation: Identifying the Euclidean function e^^ with the 

(Euclidean version of) a left moving mode, and e~^^ with a right moving 

(1) 'j'(^) 
mode, we see that r]j = -^ is the (Euclidean) reflection coefficient asso- 

dated to the left mirror, for incoming right modes. In the same fashion, 

(2) '2^(2) 

r\^ = — ^ corresponds to the reflection coefficient of the right mirror, for 

incoming right modes (for non-symmetric potentials one may have r^ ^ ^ll- 
In terms of these coefficients the Casimir force reads 

I [ d^k d ^ 



^a)^(2)^-2n. 



(31) 



which is of the same form that the Lifshitz formula for a scalar field. As 
already mentioned in Section 2, the case of finite temperature is obtained by 
replacing the integral over /cg by a sum over dicrete Matsubara frequencies. 

3 Examples 

The matrices T^*^ appearing in the final expression of the force can be explic- 
itly computed in several particular cases. For instance, let us assume that a 

10 



mirror is described by a square potential barrier of height Vo{k) and width 6 
|13] . A straightforward calculation yields: 

Til = cosh.{Q6) H —^ — sinh(r25) 

T22 = cosh.{Q6) sinh(fi5) 

T12 = -T21 = ^^~-^ smh{n6) , (32) 

where Q'^ = Q"^ + Voik). Note that for these symmetric potentials, the left 
and right reflection coefficients coincide, that is rl = r^ . 

The Casimir force between two mirrors described by such potentials is ob- 
tained by replacing the corresponding matrix elements into Eq. (l30|) . and the 
final result is in agreement with Ref.jlS]. We show explicitly the expression 
of the force only in some limiting cases: for very thick slabs with 6i » I, 
the vacuum force reads 



I r d^k d , r. V'^'Vr 



^=-7; 77^7^7 log 



(1)t7(2) 



"^0 -2f7« 



(33) 



2j {2'kY dl ^V (fi + fi(i))2(fi + l](2)) 

On the other hand, one can also obtain from the general expression (|32|) the 
T matrix associated to a singular potential of the form V{x, k) = X{k)6{x), 
where 6{x) is the Dirac 5- function. Indeed, taking the limits Vo{k) ^ 00 and 
(5—7-0 with Vo{k)6 — t- X{k), the T matrix becomes: 

T = ( ^^f 1 ^A ) ' ^^^^ 

thus, the force between two zero- width mirrors, with functions Xi{k) and 
X2{k), at a distance / appart is 



1 f d^k d ^ ^ ^^ 



(2^)2 ^_2ni 



^ ^ 20^"^-^ ^ 20' 



(35) 



which reproduces a well known result [T0| [T3| ITi] . 

4 Conclusions 

In this paper we presented a simple derivation of Lifshitz formula, based on 
the use of GY theorem to compute the functional determinant of an operator 
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of the form —d"^ + rn^ + V , where \^ is a background potential that models 
the imperfect mirrors. 

The fact that the Casimir force between flat mirrors can be written just 
in terms of the reflection coefficients of each mirror, is quite transparent in 
this derivation. Indeed, after separation of variables the problem reduces to 
the calculation of the functional determinant of the one dimensional operator 
T{k) = —d'j + Q'^{k) + V{xd, k). GY formula implies that this determinant is 
proportional to ip{L/2), where ip{xd) is a solution of T{k)ip{xd) = 0, satisfy- 
ing the initial conditions ■?/;(— L/2) = 0, ip'{—L/2) = 1. The key observation 
is that one can evaluate ip and its derivative on the right side of each mir- 
ror in terms of their values on the left side, by means of a transfer matrix 
determined by the potential describing the mirror. This matrix contains all 
the information of the potential upon which the functional determinant may 
depend. On top of that, the Casimir force only depends on some particular 
ratio of elements of the transfer matrix. Hence, the apparent dependence on 
the function that describes the mirror collapses to a single number (which 
may of course depend on frequency or momentum). 

We wish to point out that the self-energies of the mirrors do not appear in 
the expression for the force, since the derivative with respect to the distance 
between mirrors erases them out. The same would happen if we considered 
the difference between the energies corresponding to two different distances. 

Finally, we note that expression fl30|) for the force can be thought of 
as a generalization of Lifshitz formula to the case of mirrors described by 
potentials of the form V = V{xd, Xq — x'q, x\\ —x'n), that may include not only 
potential barriers but also potential wells. In this case, equation (130|) holds 
true, even when its interpretation in terms of reflection coefficients in the real 
time formulation would be problematic, because of the possible existence of 
bound states. Moreover, potentials of this form are generated in a wide class 
of theoretical models in which the vacuum field is coupled to the internal 
degrees of freedom of the mirrors, and therefore this approach provides an 
explicit link between the force and the microscopic interactions. 
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